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ABSTRACT 

Transiently accreting neutron stars in quiescence (L^^IO^'' ergs""'^) have been ob- 
served to vary in intensity by factors of few, over timescales of days to years. If the 
quiescent luminosity is powered by a hot NS core, the core cooling timescale is much 
longer than the recurrence time, and cannot explain the observed, more rapid variabil- 
ity. However, the non-equilibrium reactions which occur in the crust during outbursts 
deposit energy in iso-density shells, from which the thermal diffusion timescale to the 
photosphere is days to years. The predicted magnitude of variability is too low to 
explain the observed variability unless - as is widely believed - the neutrons beyond 
the neutron-drip density are superfluid. Even then, variability due to this mechanism 
in models with standard core neutrino cooling processes is less than 50 per cent - still 
too low to explain the reported variability. However, models with rapid core neutrino 
cooling can produce variability by a factor as great as 20, on timescales of days to 
years following an outburst. Thus, the factors of ~few intensity variability observed 
from transiently accreting neutron stars can be accounted for by this mechanism only 
if rapid core cooling processes are active. 
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1 INTRODUCTION 

X 
S3 



accretion in quiescence (Van Paradijs et al. 1987), possibly 



Many i eutron star (NS) X-ray binaries go tlirough accre- 



tion outbursts (_Lx~lU ergs 



tollowed by long periods 

,34 „„„-l\ 



(months-decades) of relative quiescence (Lx^lO" ergs 
The origin of these outbursts remains under debate, al- 
though many agree that in the wider binaries an accretion 



through an advection-dominated accretion flow (Narayan 
et al. 1997; Menou et al. 199E); accr etion onto the NS 
magnetosphere (Campana et al. 1997) following a "pro- 
peller phase", in which the NS magnetosphere is larger 
than the keplerian orbit with an orbital period equal to the 



disk in itability is the cause ( Van Paradijs 199(:). During 
these outbursts, on the order or^Tn''''g'^nay be accreted, 
over a pe riod of --^few-SO days (for a review of these tran- 
sients see Tanaka fc Shibazaki 199( ; Campana et al. 1997 



spin period of the NS (lllarionov & Sunyaev 1975 



Stella 



Chen et al. 19971) 



Following an outburst, the NSs return to quiescence, 
and have typically been detected with luminosities in the 



1032_io33 ergs"^ range (Van Paradijs et al. 1987 


IVerbunt 


et al. 1994; 


Asai et al. 1996a; Asai et al. 19961; 


Campana 


et al. 1998a 


Garcia & Callanan 1999; Rutlcdge et al. 19991 


Campana et al. 2000 


Rutledge et al. 2000). The quiescent X- 



etaL198£); and to thermal emission from a hot NS core, 
heated by non-equ ilibrium reactions deep in the NS crust 
(Brown et al. 1998, BBR98 hereafter). Of these possibilities, 
only deep-crustal heating predicts the similar luminosities 
observed in the thermal spectral component of these sys- 
tems. The observed temperatures of the thermal component 
(fcT=0.08-0.20 keV) are as expected in this scenario. In ad- 
dition, the emission area radii of the thermal component are 
consistent with theoretically predicted NS radii (~10 km; 



Rutledge et al. 1999; Rutledge et al. 2000); this supports 



ray spectrum has been described as being comprised of two 
components, a soft (black-body kT ~0.2 keV) thermal com- 
ponent, and a po wer-law component dominating the emis- 
sion above 2 keV (|Asai et al. 1996b Campana et al. 1998b 



Campana et al. 200C ) 



The basal luminosity has been attributed to continued 



the interpretation of this spectral component as a thermal 
NS photosphere. It is possible that the quiescent luminosity 
is due to a combination of these three mechanisms, in which 
case they are only separable by their spectral and intensity 
variability properties. 

The detected flux during quiescence from some of these 
transients has been observed to vary, by a factor of as much 
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as ~3 or more, over timescales of days to years. No clear 
cause of this variability has yet emerged. There are three 
means discussed in the literature: 



(i) Variable accretion onto the NS surface. If accretion 
onto the NS surface dominates the luminosity in quiescence, 
changes in the accretion rate could account for variations 
in the thermal (fcTeff~0.2 keV) part of the spectrum, since 
the photospheric spectrum at the low implied accretion rat e 
is thermal in the absence of a shock (Zampieri et al. 1995), 
assuming quasi-steady-state accretion at the appropriate ac- 
cretion rates. 

(ii) Variable accretion onto the NS magnetosphere. The 
power-law spectral "tail" which can dominate emission 
above 2 keV may be due to a ccretion onto a NS magneto- 
sphere (Campana et al. 1998a). Variations in accretion rate 
onto the magnetosphere would cause intensity variability in 
this power-law component. 

(iii) Variable absorption column density (Nn)- Outbursts 
may be accompanied by outflows, increasing locally the col- 
umn density of absorbing material. This would cause vari- 
ations in the absorption column density (A''h), which may 
account for the factor s of 2-3 variation in intensity observed 
(iRutledge et al. 199£|). 



Moreover, intensity variability due to accretion - either pho- 
tospheric or magnetospheric - would likely be stochastic; 
there have been no theoretical estimations of its magni- 
tude or variability timescale. In the variable- A''h scenario, 
the amount of absorption should decrease (and intensity, in- 
crease) with time after outburst. 

In this paper, we investigate the time dependence of the 
thermal emission mechanism proposed by BBR98. Compres- 
sion of the NS crust by accretion during outbursts induces 
electron captures, neutron emissions, and pycnonuclear re- 
actions deep in the crust (at densities > 10^ g cm~'^), and 
these reactions deposi t Qnuc ~ 1 MeV per acc reted baryon 
of heat into the crust ([Haensel fc Zdunik 1990h|). BBR98 ar- 



gued that these reactions heat the NS core to an equilibrium 
temperature ~ 10* K, and, during quiescence, the hot core 
shines with a typical luminosity Lq « {M){Qnuc/'mi,), where 
(A/) is the mean accretion rate, averaged over the thermal 
time of the core, i.e., many recurrence intervals Tree. Con- 
sequently, the core luminosity is not expected to change on 
timescales shorter than its thermal time, i.e., ~ 10^ yrs. 

Transient X-ray emission from the non-equilibrium re- 
actions has not been previously considered. Similar work 
regarding transient energy deposition in NSs has either fo- 
cussed on the atmosphere, well above the crust ( [Eichler "& 



et al. 1997 



Cheng 1989) . or investigated the thermal response to pul 
sar glitches ( vaji^jli£e£j_991 ; Chong & Cheng 1994; 



Hiranc 



Cheng et al. 1998). However, the crucial differ- 



ence is that the depths and amounts of energy deposited 
by non- equilibrium reactions can, at least in principle, be 
obtained from an ab initio calculation (and have been, for 



the case of iron, by Haensel & Zdunik 1990b), while, on the 



other hand, the total amount of energy deposited by the 
crust-breaking glitches is not precisely known (it is bounded 
by 210,60,), and neither the depth nor the distribution of 
energy deposition is well constrained. 

Colpi et al. (2001) examined the effect of non- 



~1 year, the core heats to an equilibrium temperature 
in ~10^ yrs, after which the core temperature varies by 
~0.5% in response to the periodic input of energy from the 
non-equilibrium reactions during individual accretion out- 
bursts. They confirmed the steady-state assumption made 
by BBR98. Our study is complimentary to that of Colpi 
et al., as we examine the time dependence of the tempera- 
ture in the crust and the thermal luminosity from the NS 
surface between outbursts. We follow in detail the thermal 
relaxation of the crust and find that, depending on the mi- 
crophysics and the accretion history, the magnitude of the 
variability can be as small as < 1%, or as large as a factor of 
20. Only the reactions that deposit energy at depths where 
rth<Trcc lead to variable thermal emission. The luminosity 
due to individual reactions is largely blended together, sim- 
ply because the difference between the Tth of adjacent reac- 
tions is typically much shorter than rth- However, the emis- 
sion from the reactions in the outer crust as a whole is typ- 
ically well-separated from that due to the reactions in the 
inner crust, resulting in characteristic "double- hump" lumi- 
nosity evolution. In addition, the amplitude of the variability 
is not simply related to the energy deposited in the individ- 
ual reactions. We quantify the dependence of this variability 
on the conductivity of the crust and the physics of neutrino 
emission in the core. 

In §^ we briefly recount the observational evidence for 
quiescent intensity variability. In §^ we describe detailed 
calculations of time evolution in the thermal NS flux, in re- 
sponse to a "delta-function" (1-day long) accretion event, for 
different time-averaged accretion rates (and thus core tem- 
peratures), different dominant conductivies and core neu- 
trino cooling prescriptions We also present the resulting light 
curves in this section. Finally, in we compare these re- 
sults with observations, and in 5H, we conclude. 



2 OBSERVATIONS OF VARIABILITY IN 
QUIESCENCE 

Intensity variability in quiescent transients is, at present, not 
well-studied observationally, largely due to the low signal-to- 
noise of the data, from a small number of observations (only 
a few of which are made with the same instrumentation), 
and different assumed intrinsic energy spectra. Multiple ob- 
servations in quiescence (at Lq<^W^'^ ergs~^) are found in 
the literature only for three NSs, Cen X-4, 4U 21294-37, 
and Aql X-1. 



equilibrium reactions on the temperature of the core. They 
found that in transient NSs with recurrence timescales of 



For Cen X-4, Van Paradijs et al. (1987) found the lumi- 
nosity increased by a factor of ~ 2—5 over ~5.5 yr (Jul 1980- 
Feb 1986), for the same assumed thermal bremmstrahlung 
spectrum. There was no intervening outburst observed, and 
the observations were made using two different instru ment s 
{EXOSAT/LE and Einstein /IPC). Campana et al. ( |l997| ) 
reanalysed these data, and contrarily concluded that they 
are consistent with the same luminosity. Also, Campana 
found that in observations with ROSAT/RRl of Cen X-4 
over a 4-8 day period, the source countrate varied by ~3, 
with an average luminosity of ~7xl0^^ ergs~^; Campana 
did not discuss whether the 1995 ROSAT/URl observa- 
tions were consistent with the luminosities of the EXOSAT 
and Einstem observations; a comparison between their val- 
ues and passbands indicates that the luminosities of Van 
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Paradijs et al. were more than a factor 3 greater, although 
part of this may be du e to a different assumed spectrum. 
Finally, Rutledge et al. (2001) found that Cen X-4 varied by 
<18% (0.2-2 keV) during a 10 ksec observation, while the 
0.5-10.0 keV luminosity had decreased by 40±8% in com- 
parison with an ASCA observation taken 5 years previously. 

Aql X-1 was also found to have a variable flux in qui- 
escence (Rutledge et al. 1999) between three observations 
with ASCA (Oct 1992, Mar 1993, and Oct 1996). If the 
same intrinsic spectrum is assumed, the three observations 
are inconsistent with having the same intensity with high 
confidence (probability=10~®). The difference in observed 
flux is a factor of ~2-3. 

4U 2129-1-37 decreased in flux by a factor of 3.4±0.6 
between Nov-Dec 1992 (ROSAT/URl) and M a rch 1994 
(ROSAT/PSPC) dCarcia fc Callanan 1999|; [Rutledge 



et al. 2 )0C ) in the unabsorbed luminosity, for the same as- 
sumed spectrum. 

The observations of Cen X-4, with ROSAT/URl, and 
of Aql X-1, with ASCA are the only two instances of re- 
peated observations with the same instrumentation. C en X- 
4 has n ot b een observed in outburst since May, 1979 (Chen 
et al. 1397), while Aql X-1 goes into outburst every ~220 
days, and it therefore seems likely that there were interven- 
ing outbursts between the three ASCA observations. 

Thus, observations of transiently accreting NSs in qui- 
escence, taken at face value, indicate that their quiescent 
luminosities vary by a factor of up to 3—5 on timescales 
of days to years. However, it is not clear what fraction of 
this variability is intrinsic, and what can be attributed to 
systematic differences in instrumentation. 



3 THERMAL EMISSION FROM THE CRUST 

We now describe our simulations of the thermal relaxation 
of the crust following an accretion outburst, and the result- 
ing time dependence of the quiescent thermal e missio n. As 
argued by BBR98 and confirmed by Colpi et al. ( 2001 ), the 
nuclear energy release in the crust heats the NS core to a 
temperature corresponding to steady accretion at the corre- 
sponding time-averaged accretion rate. We performed sev- 
eral simulations starting with a very cold NS and subjecting 
it to a series of accretion outbursts and found that indeed, 
the core is heated to the appropriate tempera ture o n a 10* yr 
timescale, in agreement with Colpi et al. (2001). For the 
remainder of our analysis we first construct a steady-state 
thermal model corresponding to the time-averaged (M), and 
then subject it to outbursts. After several outbursts, the 
model reaches a limit cycle. In this paper, we report the re- 
sults of only these latter simulations. Our accretion events 
("outbursts") are 1-day long "delta- functions" , with no ac- 
cretion outside of these events. More realistic time-evolved 
light-curves may in principle be found by convolving our 
"delta-function" thermal response light curves with a more 
realistic outburst accretion time-dependent profile for the 
small perturbation light curves [SL / L<X); lightcurves in 
which the variations are larger would require a more de- 
tailed calculation. We simulate two different outburst recur- 
rence timescales (1 yr and 30 yr) and three time-averaged 



surface flux from the crust and core to a "delta-function" 
accretion profile, it is unimportant that the implied 1-day 
accretion rates may be super-Eddington, and we do not in- 
clude accretion luminosity in our results. 



3.1 Microphysics of the Crust and Core 

The hydrostatic and steady-state thermal models of the 
crust used in this paper are substantially ident ical t o the 
ones used by Ushomirsky, Cutler, & Bildsten (2000). We 
briefiy recount the major ingredients of these models, and 
then describe our time-evolution code. 

Since the thermal timescale for the crust (days to years, 
see below) is so much longer than the sound crossing time 
(milliseconds), and the equation of state is nearly indepen- 
dent of temperature, the hydrostatic structure of the crust is 
effectively decoupled from its thermal evolution. Therefore, 
one can solve the thermal evolution equations assuming that 
the density and pressure as a function of the position in the 
star are fixed. This approach to modeling the thermal state 
of a neutron star is fully described by Brown ( 2000| ) and is 
followed here. 

The composition of the crust (i.e., the mass A and 
charge Z of the nuclei, as well as the neutron fracti on Xn) 



is take n from the tabulation of Haensel & Zdunik (1990b 



1990a). The pressure is the sum of contributions from de- 



generate, relativ istic electrons and free neutrons (Negele & 
Vautherin 197J) at de nsities exceeding neutron drip (pnd ~ 
(4 - 6) X 10^^ g cm"^, iHaensel & Zdunik 1990a|). With the 



equation of state as described above, we solve the Newtonian 
equations of hydrostatic balance (which, for a thin crust, 
are exactly equivalent to the fully relativistic equations). 
We take M = 1.4M0, and 7i = 10 km. Our crust has a 
thickness of 1.1 km and a mass of O.OGMq. 

With the hydrostatic structure specified, we solve the 
heat equation, 



^ ^ dt r'^ dr \ dr 



+ p(enuc — Ei/), 



(1) 



where c„ is the heat capacity per gram (which has contri- 
butions from the ionic lattice, degenerate electrons, and free 
neutrons in the inner crust), e^uc is the nuclear energy re- 
lease, is the neutrino emissivity, and K is the thermal 
conductivity. We neglect the overall downward motion of the 
material due to accretion, as well as the heat release due to 
compression (these effects are not important in the crust). 
Therefore, the only energy source in the crust is the heat 
release due to non-equilibrium electron captures, neutron 
emissions, and pycnonuclear reactions. The detailed treat- 
ment of the energy release is described in Ushomirsky et al. 
( ^000| ); we note here that enuc is proportional to the instan- 
taneous mass accretion rate M. The total energy released in 
the crust is taken to be Qnuc = f A-kt^ pe^i y^cdr = 1.45 MeV 



per accreted baryon ( Haensel & Zdunik 1990b ) . At the tem- 



accretion rates (10 



10' 



and 10 Mq/jv). As we are 



primarily interested in illustrating the response of the NS 



peratures of interest, neutrino losses in the crust are negli- 
gible. 

Instead of fully modeling the thermal evolution of the 
NS core, we presume that it is isothermal, and is char- 
acterised by a single temperature Tcore which is equal to 
the temperature of the crust at the interface with the core, 
T{r core)- The core temperature obeys 
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Cci 



rir(rco..c 
dt 



dT 
dr 



(2) 



where Ccoro is the total heat capacity, and L^, is the total 
neutrino emissivity of the core (see below). This approxi- 
mate treatment is justified because of the core's large ther- 
mal conductivity and heat content. In effect, the tempera- 
ture of the core changes by a negligible amount compared 
to the temperature variations in the crust, and the core acts 
as an energy sink at the bottom edge of the crust. With 
this approximation, we do not need to integrate the heat 
equation in the core, and instead use Eg. (|2[) as a boundary 



condition for Eq. (1) at the bottom edge of the crust. 



To survey the parameter space, we use two difi'erent 
prescriptions for cor e cooling. "Standard" co o ling models 
use mo dified U rea ( Friman & Maxwell 197!: ; [Vakovlev & 



Levenfi ih 1995 ) and e-e neutrino bremsstrahlung (Kaminker 
& Haer sel 199E| ), while "rapid " cooling models assume t he 
presence of a pion condensate (Shapiro & Teukolsky 1983 )Q 



Nucleon superfluidi ty significantly reduce s neutrino emissiv 

for an in-depth re- 



Yakovlev et al. 1999 



ity in the core 

view). In our calculation we use super fluid parameters col- 



lated by Brown (gOOC|). Maxwell ( |l979[ ) argued that at tem- 
peratures well below the superfluid transition temperature 
Tc, this suppression is exponential. Howe ver, d etailed cal- 



culations, summarized by Yakovlev et al. (1999) show that 



this suppression is weaker by a power-law factor in T/Tc, 
leading to lower core temperatures than one would expect 
based on the simple exponential reduction in emissivity. 

Core heat capacity is the sum of contributions from elec- 
trons, neutrons, and protons. However , superfluidity alters 
the heat capacity of the nucleons (see Yakovlev et al. 1999 
for an in-depth review). At T <^ Tc, the reduction of the 
heat capacity is simila r to the suppression of neutrino emis- 
sivity. Maxwell (1979) gives a fitting formula for ^So super- 



fiuid (applicable to neutrons in NS crust and protons in the 
core), while Levenfish & Yakovlev (1994) give fitting formu- 
lae for the superfiuid as well (applicable to neutrons 
in the core). At temperatures of interest 10^ — 10** K) 
electrons dominate the heat capacity. 

Throughout the crust, heat is transported by degener- 
ate, relativistic electrons, and the thermal conductivity K is 
very sensitive to the purity of the cr ust. If the crust is a pur e 
crystal, electron- phonon scattering (Baiko & Yakovlev 199E) 



sets the conductivity. However, crusts of accreting neutron 
stars are unlikely to be pure crystals, since they are com- 
posed of (likely quite impure) products of nuclear burning 
in the uppe r atmosphere. In fact, calculations of Schatz 
et al. (1999) suggest that the typical values of the im- 
purity parameter Qimp = Yimp{Zimp - (Z))^, where Zi^p 
and Yimp are the charge and the fraction of the impurities, 
could be comparable to the average charge (Z)'^. Follow- 
ing Brown (2000), we set the absolute lower b ound on the 



conduc tivity by using electron-ion scattering (Yakovlev & 
Urpin ] |980| ), which has the same form as electron- impurity 
scattering with Qimp = Z'^. The much smaller conductivity 
in the case of electron-impurity scattering with large Qimp 
leads to much longer crustal thermal time, and much smaller 



* Instead, one could assume that the NS mass is high en ough to 
allow direct Urea neutrino emission (see, e.g., Colpi et al. (2001)), 
as the emissivities are similar. 



variability of quiescent thermal emission. This approach al- 
lows us to survey the entire range of possibilities for the 
crustal conductivityjT] 

At the top of the crust (p ~ 10^ g cm^'^) we utilize the 



results of Potekhin et al. (1997) for the relation between 
the core temperature to the surface temperature.^ They 
extended the original calculations of Gudmunsson et al. 
( 1983| ) by considering the opacities of accreted envelopes, 
rather than iron envelopes. Because of the smaller opac- 
ity of a light-element envelope, a ~50% smaller core tem- 
perature is necessary to carry a given fiux (see also 



Bland- 



ford et al. 1983). The thermal time at the top of the crust. 



where the outer boundary condtion is applied, is 1 day in 
our models, so our simulations cannot follow variability on 
shorter timescales. 



3.2 Time- Averaged Quiescent Luminsity 

Before describing the results of our time evolution simula- 
tions, we discuss the thermal luminosity of neutron stars 
assuming steady accretion at the rate (M) corresponding 
to the time average over the outburst recurrence interval 
Tree. This provides a reasonable estimate of the average lu- 
minosity level of neutron star tra nsients. The variations of 
the luminosity discussed in § |3.3| are excursions about this 
average level. As outlined in § |3.l| , we survey the parameter 
space by considering two different strengths of core neu- 
trino emission (standard and enhanced), and two different 
assumptions regarding the conductivity of the crust. 

The average core temperature and luminosity of the NS 
are set by the balance between the heat input due to non- 
equilibrium reactions in the crust (at the rate Qnuc{M) /mb) 
and heat loss due to neutrino emission from the core (L^/) 
and photon luminosity from the surface (L-,). At the low core 
temperatures characteristic of NS transients (Tcoro <^ 10* K) 
modified Urea neutrino emission (standard cooling case) is 
substantially suppressed by nucleon superfiuidity and can- 
not compete with photon losses from the surface of the star, 
Li, <^ L~f. Therefore, most of the heat deposited by the 
nonequilibrium reactions is radiated from the surface, i.e., 
the quiescent luminosity is just (BBR98) 



(3) 



= 8.7 X lO-" 



(M) 



10-10 Mq yr- 



1.45 MeV 



erg s 



This estimate is independent of the conductivity of the 
crust. The core temperature, however, depends on the 
crustal conductivity. If the crustal conductivity is large 
(i.e., set by electron-phonon scattering), the crust is nearly 

t The presence of impurities does not dramatically affect the heat 
capacity of the ionic lattice. While the phonon spectrum will be 
affected by the presence of impurities, the heat capacity should 
still approach the Debye limit. 

t Potekhin et al. ( ) express the surface temperature in terms 



of the temperature in the isothermal region of the crust, which 
they assume to be at p > lO^^ g cm" ^. However, since the au- 



thors assume that the crust has high thermal conductivity, the 
temperature difference between lO^" g cm"^ and 10* g cm~'^ is 
negligible in their models. 
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Figure 1. Time evolution of the temperature of the crust for a model with standard cooling and low thermal conductivity crust following 
a 1-day outburst and (A/) = 10~^''Mq yr"'^, with a recurrence time of 30 years. The x axis is the density, covering the entire crust of 
the neutron star; y axis is the time, and z axis is the temperature. 



isothermal, and the core temperature, Tcorc ~ 1-2 x 
10** K ((M>/1O-^"M0 yr-'^f-^\ can be obta ined by sim- 
ply using the relation of Potekhin et al. (1997). If the crust 
has low thermal conductivity, then a substantial tempera- 
ture gradient is needed to carry the heat from the inner 
crust to the surface, and the core temperature is Tcoro ~ 
2.2 X 10** K ((A/)/lO-^°M0 yr-i)° '*'5 (obtained by fitting 
our detailed calculations). 

On the other hand, if rapid cooling processes are al- 
lowed in the core, then, despite the suppression of neu- 
trino emission by superfiuidity, most of the heat deposited 
by non-equilibrium reactions is radiated away by neutri- 
nos, Li, ^ Ly. The core temperature is then set by bal- 
ancing the energy input from the crustal reactions with 
neutrino-luminosity. Except at very low accretion rates (<; 
5 X lO~^^Af0 yr~^), where superfluid suppression of neu- 



trino emission makes comparable to L^, a good fit 
to the results of our calculations for the particular case 
of enhanced neutrino emission from a pion condensate is 
Tcoro ~ 2.7 X 10'' K ({M)/1O-^°M0 yr-^)°■°^ regardless of 
the conductivity of the crust. The thermal luminosity, how- 
ever, depends on the relation between the core temperature 
and that at the top of the crust. When the conductivity is 
determined by electron- phonon scattering (i.e., K is large), 
the crust is ne arly is othermal, and, using the relation of 
Potekhin et al. ( |l997[ ), we find 

0.22 

erg 



(M) 



10-"' Mq yr- 



Qnuc/mt 

1.45 MeV 



which agrees well with our detailed calculations. On the 
other hand, if the crustal conductivity is low, the relation 
between the core and the surface temperature is nontrivial. 
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since an appreciable temperature gradient between neutron 
drip and the core is needed to carry the flux from the non- 
equilibrium reactions into the core. In this case, the crust 
is warmer than in the electron-phonon scattering case, and 
the thermal luminosity is 



10-10 Mo yr- 



Qnuc/mi, 

1.45 MeV 



erg s 



"(5) 



(which we obtained by fitting the results of our detailed cal- 
culations). Comparing equation with ^ and (^), we see 
that when enhanced neutrino cooling processes are allowed 
in the NS core, only less than 10% of the non-equilibrium 
energy release is radiated from the surface, while the remain- 
ing > 90% of the heat is emitted in neutrinos. Regardless of 
the crustal conductivity, the core temperature in the rapid 
cooling case is set entirely by the core neutrino emission. 
However, depending on the conductivity of the crust, a tem- 
perature gradient between the core and the surface may or 
may not be present, resulting in different surface tempera- 
tures and photon luminosities (cf. Eqs. and (^)). Thus, 
neutron stars with rapid cooling processes active in the core 
appear dimmer than those with just standard cooling, and, 
when their cores are very cold due to rapid neutrino cooling, 
NSs with high conductivity crusts will appear dimmer than 
those with lower conductivity crusts. 



3.3 Time- Variable Emission 

As we now describe, the approximate steady-state luminos- 
ity is only a part of the story. Depending on the microphysics 
of the crust and core, the time-dependent thermal luminos- 
ity of the neutron star may either always be very close to the 
steady-state estimate {SL/L < 1%), or vary wildly around 
it (SL/L^l). 

In Figure |l| we display the time evolution of the temper- 
ature in the NS crust for a model with low thermal conduc- 
tivity, standard neutrino emission, (A/) — 10~^^Mq yr~^, 
and rrc c=30 vrs. Other models are similar qualitatively, but, 



of cour ie, differ in the magnitude of variability. The slices 
of the figure in the x — z plane represent the temperature 
in the crust as a function of density, and the y axis is the 
time (in days) since the beginning of a 1-day outburst. At 
the early times, the temperature rises locally at the densi- 
ties where the energy is deposited. It reaches the maximum 
after 1 day (i.e., at the end of the outburst). In the en- 
suing cooling period, the heights of the temperature peaks 
decrease, and their widths increase due to heat diffusion. 
After ~ 6 days the "heat wave" due to the very first reac- 
tion, at p = 1.5 X 10^ g cm~^ reaches the top of the crust. 
This temperature increase at the top of the crust is then 
directly transla ted into the increase in the thermal luminos- 
ity of the NS ( |Potekhin et al. 1997| ). At much later time 
(~ 1000 days after the outburst) the heat wave due to the 
reactions in the inner crust, at p ^ 10^^ g cm~'^ reaches the 
surface. Even though the energy release in these reactions 
is much larger, the change in the outer boundary tempera- 
ture and, hence, the variation of the thermal luminosity is 
much smaller, because the temperature peak spreads out to 
a width corresponding to roughly the thermal time at the 
deposition depth as it travels towards the surface of the star, 
and because some of the heat is conducted into the core of 



the star. Moreover, while the temperature peaks in the crust 
due to individual reactions are well-separated at the early 
times, they become blended together by the time they reach 
the surface. 

Figure |^ clearly suggests that, in order to produce ap- 
preciable variability in the thermal emission from the sur- 
face, the energy deposition by the non-equilibrium reactions 
should be sufficient to heat the crust locally by ST ~ T. 
The typical outburst fluences indicate that only AAI <. 
10^^ grams of material is accreted during an outburst, and 
hence < 10*^ erg of energy is deposited in the crust. This 
amount is rather small, which allows us to place some in- 
teresting constraints on the heat capacity of NS crust. The 
thermal time at the depths where this energy is deposited 
is at least comparable to and usually much greater than 
a typical outburst duration. Thus, during the outburst, the 
heat is not conducted away, but primarily heats the crust lo- 
cally. Neglecting heat diffusion, the nuclear energy deposited 
in the region near neutron drip will heat this region by an 
amount 



ST' 



C 



ks /baryon 



Pd 



Qnuc \ f AM 



(6) 



U MeVy 



\^1023 



where C is the heat capactiy in units of ks (Boltzmann con- 
stant) per baryon and pd is the pre ssure at which the energy 
is deposited. As outlined in § |3.l| , the typical crustal tem- 
perature is ~ lO'^ — 10* K. Therefore, if the heat capacity of 
the neutron star crust is ~ fcs per baryon, the crust relax- 
ation luminosity is only a small perturbation on the overall 
cooling of the core. In order to produce observable variabil- 
ity of quiescent thermal emission, the NS crust must have 
C <^ ks per baryon. 

Figure |^ shows contributions to the heat capacity in the 
crust, for a representative case of (M) — 10~^^Mq yr~^. 
In the outer crust, the heat capacity of the ionic latice , 
C < lO'^fcfl per baryon (dotted line, Shapiro fc Teukol 



sky 1983| ; |van Riper 199ll ; |Chong fc Cheng 1994|) , dominates 



Strongly degenerate, relativistic electrons (dashed line) con- 
tribute a negligible amount to the heat capacity everywhere 
except in the deep inner crust (p ^ 2 x lO'^^ g cm^). At den- 
sities above neutron drip, free neutrons, if not superfiuid, 
would have appreciable heat capacity, approaching ks per 
baryon (dot-dashed line). Thus, if the neutrons in the inner 
crust are not superfluid, their heat capacity is so large that 
the energy release of < lO**^ erg characteristic of the heating 
during an outburst is not sufficient to heat the crust substan- 
tially. However, it is commonly believed that the neutrons in 
the inner crust for m a^Sp superfluid, and their heat cap acity 



is greatly reduced (Maxwell 1979; Yakovlev et al. 1999). Su- 



perfluidity completely suppresses free neutron heat capacity 
at densities > 8 x 10^^ g cm'' (solid line in Figure ^). In this 
case, the ions dominate the heat capacity in the inner crust, 
and C < 10~^fc_B per baryon. Therefore, observing the evo- 
lution of quiescent luminosity provides a test of superfluidity 
of neutrons: if late-time evolution is observed {ST /T >10~''), 
then the neutrons must be superfiuid. While this may not be 
a controversial conclusion, it is worth noting that the only 
existing observational evidence for superfiuidity of crustal 
neutron gas is the interpretation of pulsar glitches as be- 
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Figure 2. Heat capacity of neutron star matter as function of 
depth, for (M) = 10~^^ Mq yr~^, in units of kg per baryon (in 
these units, the heat capacity of ideal monoatomic gas would be 
C = 3/2A). Dashed line: degenerate relativistic electrons, dot- 
ted line: ionic lattice contribution; solid line: degenerate, non- 
relativistic, superfluid neutrons; dash-dotted line: degenerate non- 
relativistic normal neutrons. Note that the heat capacity of de- 
generate neutrons first rises, and then drops dramatically at the 
onset of superfluidity, at p ^ 7 X lO^'^g cm~^. The heat-capacities 
at (M) =10"-'^'' and 10"-*^^ Mq yr"-"^ are qualitatively similar. 



ing initiated in the crust. Variability of quiescent thermal 
emission from NS transients, if disentangled from accretion 
in quiescence or other sources of variability outlined in § |l[ 
would be a completely separate direct observational confir- 
mation of this view. 

The spreading and blending of the temperature peaks 
due to different reactions, clearly evident in Figure hJ, as 
well as the relative magnitude of heat diffusion towards the 
surface and the core can be easily understood in terms of 
the thermal time at the place where the heat is deposited. 
To make the discus sion more precise, we defin e the thermal 
diffusion timescale ([Henyey fc L'Ecuyer 19691; BBR98) 



Tth = 



1/2 



dr 



(7) 



Roughly speaking, rth is the time it takes a heat impulse 
to diffuse from ri to r2. In Figure ^ we plot the thermal 
diffusion timescales as a function of density in the crust to 
the NS surface (solid line) and the core (dotted lines). The 
top panel is for the model with electron-impurity scattering 
conductivity in the limit of very impure crust, while the bot- 
tom panel is for the model with electron-phonon conductiv- 



S 

ID 



03 

B 




Density (g cm ^) 



Figure 3. Thermal diffusion time to the surface (solid lines) and 
to the core (dotted lines) as function of density for (reading lines 
top to bottom) (A/) = 10-10, 10"", and IQ-^^ Mq yr"!. The 
top panel corresponds to models with standard core cooling and 
low crustal thermal conductivity, while the bottom panel is for 
models with standard core cooling and high crustal thermal con- 
ductivity. The filled bars in the top panel indicate schematically 
the regions of nuclear energy deposition. The height of the bars 
is proportional to the deposited energy, with the highest bar at 
p 1.5 X lO^'^ g cm~^ corresponding to 0.47 MeV per accreted 
baryon deposited. 



ity.^ Black vertical bars in the top panel show the locations 
of the energy deposition due to non-equilibrium reactions, 
with the height of the bars proportional to the energy de- 
posited. Clearly, the energy depositions in the inner crust 
happen at very similar depths, such that the differences be- 
tween the thermal times of the reactions are much smaller 
than the thermal times themselves. By the time the heat 
from these reactions diffuses to the surface, the differences 
between them are blended. It is also clear from this figure 
that if the thermal diffusion time to the core is smaller than 
that to the surface, then most of the heat will flow to the 
core, rather than to the surface. 

In Figures ^ and ^ we survey the dependence of the 
transient emission following an accretion outburst on various 
model parameters. Figure H shows results for a crust with 



3 The thermal times are shorter for the enhanced cooling models 
(not shown here) owing to the lower crustal temperatures. 



© DODO RAS, MNRAS 000, 000-000 



8 Greg Ushomirsky and Robert E. Rutledge 



low thermal conductivity (electron-impurity scattering with 
Qimp — Z'^), while Figure ^ shows results for high thermal 
conductivity (electron-phonon scattering). These figures are 
plotted with luminosity normalized by the thermal luminos- 
ity just prior to the outburst, which is the lowest luminosity 
of the outburst cycle. This tends to exaggerate visually the 
range of luminosities. It also does not display the absolute 
differences in the luminosities due to different {M), which 
are discussed in § However, it does make clear the rel- 
evant factors of variability which can be produced over a 
recurrence cycle due to this mechanism, and that is our goal 
here. 

As is clear from the figures, the variability is the small- 
est (< 1%) for models with standard cooling and short (1 
year, in our case) outburst recurrence times. The variabil- 
ity is largest for models with rapid core cooling and long 
(30 year) outburst recurrence time. Under the same condi- 
tions, the variability of the models with high crustal ther- 
mal conductivity (Figure Bl) is always larger than that of the 
low conductivity models (Figure All of these features are 
straightforward to understand in light of the above discus- 
sion. Since the energy deposited into the crust is not large 
(< 10*^ erg), the colder the crust, the larger the transient 
thermal response. Models with rapid core cooling and high 
thermal conductivity have systematically lower crust tem- 
peratures (see § 3.^), and hence have larger amplitude of 
variability. 

As discussed above, the first hump in the light curve is 
due to reactions in the outer crust, while the second hump 
is due to those in the inner crust. This is a generic feature 
of all models except for the low-ZsT, short-recurrence ones 
(left column of Figure ^) . As is evident from the top panel 
of Figure ^, the thermal diffusion time from the inner crust 
to the surface is > 10"^ days, i.e., longer than the recurrence 
time. Thus, the second peak in the lightcurve does not make 
it to the NS surface before the next outburst, and only one 
peak is observed in the quiescence lightcurve. This mismatch 
between the thermal diffusion timescale and the outburst re- 
currence timescale explains why only one light-curve peak is 
observed in the Troc^l yr, low thermal conductivity simula- 
tions (Figures ^a & b) while two peaks are observed in all 
other simulations. For long Tree sources, the time between 
the outburst and the second peak in the lightcurve is an 
excellent observational indicator of the thermal time in the 
inner crust. 

As we see in Figure ^ the thermal time of the inner 
crust in low- if models is « 10^ — 10^ days, while it is only 
~ 10^ days for the high-Tf models. As a result, for \aw-K 
models with rroc = l yr, the crust does not completely cool 
down, and therefore remains at a roughly constant tempera- 
ture. On the other hand, for Tree =30 yr, the crust does have 
enough time to cool down, and hence the variability is much 
greater (tens of per cent to factors of few). In addition, it is 
clear from Eq. (^ that it is AM = (M)rroc, rather than {M) 
or Tree individually that is important for the overall magni- 
tude of variability. For a given (M), models with larger Tree 
exhibit larger variability, simply because the amount of en- 
ergy deposited into the crust is larger. This explains why the 
rrec = l yr light curves are (relatively) less variable than the 
Troc=30 yr light curves. In general, one therefore expects a 
greater magnitude of (relative) intensity variability in quies- 



cence from systems with longer recurrence timescales (^10 

yr)- 



4 COMPARISON WITH OBSERVATIONS 

Is the phenomenon we describe here observed in transient 
NSs in quiescence? If so, both the time-averaged luminos- 
ity level, and the variability around it need to be self- 
consistently accounted for. This constraint provides an ob- 
servational discriminant for the various core cooling scenar- 
ios in NS transients. 

The observed magnitude of variability (factors of ~ 
3 — 5) cannot be achieved in models with standard neu- 
trino cooling, which vary by less than 50% in all cases we 
investigate here. On the other hand, rapid neutrino cooling 
models are capable of producing intensity variability by fac- 
tors of up to 20 on timescales of months-years; this may be 
invoked for the observed variability of 4U 2129-1-37 (a factor 
of 3.4 over months), and the longer time-scale variability of 
Cen X— 4 (factor of 2-5 over 6 years, an additional factor of 
3 over 10 years). The factor of 2-3 variability in Aql X— 1 
over months-years can be similarly explained, again requir- 
ing rapid cooling to be operative in its core. However, to 
explain the factor of ~3 variation over a period of 4-8 days 
in the quiesc ent X -ray flux of Cen X— 4 observed by Cam- 
pana et al. (1997), one must invoke an unobserved, short 
(< 1 day) outburst, preceeding those observations by only 
~ few days, which seems unlikely, although it cannot be 
excluded due to the lack of X-ray coverage. 

However, as argued by BBR98, the inferred quiescent 
luminosity of Aql X— 1 can be naturally understood if the 
neutrino losses from its core are negligible, i.e. if only 
standard cooling processes operate (cf. our Eq. Q). Time- 
averaged Lq of Aql X— 1 is inconsistent with rapid cooling, 
regardless of the assumed crustal conductivity (cf. Eq. W 
and D). This is also the case for 4U 1608-52. For these 
sources, variability in quiescence must be due to processes 
other than thermal relaxation of the crust (outlined in § 

Of the observed recurrent NS transients, only Cen X— 4 
has time-averaged Lq much lower than given by Eq. (^. 
BBR98 conjectured that this discrepancy is due to the fact 
that only Qnuc ~ 0.1 MeV, rather than ~ 1 MeV is de- 
posited in its crust while Colpi et al. (2001) argued that 
Qnuc ~ 1 MeV, but rapid neutrino cooling is operating in 
Cen X— 4. Clearly, analysis of variability of Cen X— 4 in qui- 
escence (or other transients with long recurrence times) pro- 
vides a discriminant between the various core cooling mech- 
anisms. If enhanced neutrino cooling is responsible for the 



relatevely low Lq of Cen X— 4, then its quiescent luminosity 
must vary in a well-defined way (see panels d of Figures W 
and ^). In the contrary case of standard cooling with low 
Qnuc, the baseline thermal luminosity of Cen X— 4 should 
be constant to better than 50%, and any variability must be 
accounted for by the processes outside the NS crust. Finally, 
it is difficult to reliably predict the expected quiescent lumi- 
nosity of 4U 2129-1-47, due to its unusual accretion history 



^ The low quiescent luminosity could also potentially be ex- 
plained by an overestimation of the time-average accretion rate 
of Cen X— 4, due to the small number of outbursts observed. 
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Figure 4. Thermal emission luminosity from the NS photosphere, due to crustal heating and NS core temperature (neglecting accretion 
luminosity) as a function of time since a one-day-long "outburst" - for crustal conductivity due to electron-impurity scattering. The 
luminosity is normalized by its value just prior to the outburst (i.e., by the lowest value over the recurrence time). Panels a and c are 
for models with modified Urea ("standard") cooling, while panels b and d are for models with accelerated cooling (in this case, pion 
condensate). The assumed outburst recurrence time is indicated in each panel. The time average accretion rates used are: (M) = 10"'^'' 
(solid line), 10"^^ (dotted line), 10"^^ (broken line) Mq yi~^ . 



sions fo^ this source. 



( Pietsch et al. 1986 ), and so we draw no additional conclu- 



If transiently accreting pulsars in quiescence fully ex- 
clude accretion onto the NS surface due to the propeller 
effect, then observations of the thermal pulse during quies- 
cence will be due only to crustal and core emission. More- 
over, optical observations can indicate when the acc retion 
disk is absent in these sources (Negueruela et al. 2000), and 



observations which follow the pulsed intensity profile of such 
an object over several years following an outburst can read- 
ily constrain the NS atmospheric luminosity. The quiescent 
emission of transient pulsars, such as the pulsed emission 
recently detected from A0535-f35 (Negueruela et al. 200C) 



and the persistent emission from SAX J1808.4— 3658 (Stella 



et al. 2QQC; Dotani et al. 2000), can potentially be used for 



measuring the late-time emission predicted here. However, 
theoretical modeling of the quiescent emission in this case 
is complicated by the need to understand the effects of the 
magnetic field, transverse heat flow in the crust, and whether 
or not the crust is replaced by accretion in these systems. 



5 CONCLUSIONS 

We have investigated the time- variable luminosity of a tran- 
siently accreting NS, following a "delta- function" -like accre- 
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Figure 5. Same as Figure^ except the assumed crustal conductivity is due to electron-phonon scattering. 



tion outburst (one day in duration), for three different values 
of (M) , for NS cores with both standard and rapid core cool- 
ing (specifically, for u emissivity from modified Urea and 
from pion condensate), and for high crustal thermal con- 
ductivity dominated by electron-phonon scattering (for the 
case of a pure crystal crust) or for low conductivity due 
to electron-impurity scattering with large impurity fraction 
(for the more likely case of very impure crust expected in an 
accreting NS). 

The magnitude and time-dependence of late-time emis- 
sion due to heat deposited by non-equilibrium reactions in 
the crust depends intimately on the microphysics active in 
the crust and the core, the time-average accretion rate (and 
thus the time-average quiescent luminosity itself) and - most 
importantly - on the recurrence timescale of the transient 
system. In the case of standard core cooling, the magnitude 
of the luminosity swings (peak-to-peak) is <,!% for rrcc=l 
yr, and up to 15-40% for Troc~30 yr, independent of the crust 



conductivity and (A/) . For rapid core cooling, the magnitude 
of the luminosity swings varies between ~few percent, to as 
much as a factor of 20, and is largest for long Tree and high 
crustal thermal conductivity. 

It is possible to observe either one peak in the lightcurve 
(due to the blended emission from the several reactions in 
the outer crust) or two peaks (the second is due to blended 
emission from reactions in the inner crust). Observing the 
second peak requires a recurrence time long enough to per- 
mit the inner-crust emission to reach the surface. For low 
crustal thermal conductivity, this thermal time is ~ 10 yr, 
so the second peak is only observed in transients with longer 
recurrence times. In the electron-phonon scattering domi- 
nated crust, with thermal time of 1 yr, the second peak 
may be observed in short Tree sources. In general, only the 
energy deposited at depths where the thermal time is smaller 
than the outburst recurrence interval will produce observ- 
able variability in the quiescent thermal emission. 
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The changes in luminosity which would permit one to 
follow the time- variable emission due to deep-crustal heating 
at the ~few per cent level are 5I/X ~10^° ergs~^ (=5xlO~^® 
[d/4.0kpc]'^ ergcm~^s~^), which is detectable with Chan- 
dra, XMM and Con-X in ~10^, 5x10", and 10'* sec, respec- 
tively. Complicating such observations is the contribution of 
accretion onto the NS, which could add a stochastic vari- 
ability component to the thermal spectrum which may not 
easily be distinguished. If accretion only occurs at the mag- 
netosphere, this contributes to the power-law component, 
which can be separated spectrally. 

Our comparisons between these results and previously 
reported observations of transiently accreting NSs in quies- 
cence indicate that the observed magnitude of intensity vari- 
ability is much greater than we predict for standard neutron 
star core cooling and crust conductivity (factors of ~few for 
the former vs. < 50% for the latter). Hence, for the recurrent 
transients Aql X-1, 4U 1608-52, and SAX J1808-36, whose 
time-averaged luminosity agrees well with the predictions of 
standard neutrino cooling (BBR98; cf. our § 3.2), it seems 
likely that some other process - such as continued, variable 
accretion in quiescence - is required to explain the majority 
of the observed variability. 

In the case of Cen X— 4, however, the observed quies- 



cent luminosity is consi sten t with rapid core cooling (Colpi 
al 



et 



001 



cf. 



3.2), and hence its long-timescale 



(months to years) variability could be accounted for by 

crustal relaxation .[J] On the other hand, if Cen X— 4 ap- 
pears too dim because its time-averaged accretion rate has 
been overestimated or non-equilibrium reactions deposit an 
unusually small amount of energy into its crust (BBR98), 
then the majority of its quiescent variability must be due to 
stochastic sources external to the NS. Disentangling the true 
thermal emission from Cen X— 4 from the contamination due 
to continued accretion in quiescence (or other sources) would 
allow one to discriminate these two hypotheses. 

Finally, we note that the observations are low quality 
(S/N~few); there are only two instances of repeated obser- 
vation with the same instrument, and in those cases the 
number of observations were only 2-3 (not densely sam- 
pled light curves); and the infered time-average accretion 
rate (and, hence, the predicted quiescent luminosity) of Cen 
X-4 is highly uncertain due to the small number of out- 
bursts observed. Therefore, a definitive statment on this 
subject requires greater observational scrutiny, with higher 
S/N data than presently available, and with densely sampled 
lightcurve following an outburst, which would be capable 
of distinguishing between the temporal variability predicted 
here and more stochastic variability which likely accompa- 
nies a variable accretion rate. 
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